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1. Introduction 



Many-body condensed matter physics, as opposed to few body nuclear physics, 
molecular physics, or quantum chemistry, is concerned with the limiting case N — > oo, 
where N is the number of particles of the system. There are, of course, a large number 
of physical questions in which we are interested for these systems, such as response 
functions at finite temperature, but conceptually the easiest and most primitive ques- 
tion is the solution of the N body Schroedinger equation in the N — » oo limit. In this 
study it is surely conceptually easier to consider the spectrum of eigenvalues before 
one investigates the more detailed information contained in the eigenfunctions. Thus 
the first most primitive question asked about such systems as the free Bose or Fermi 
gases, the Landau theory of the Fermi liquid, and the BCS theory of superconductivity 
is how to characterize the eigenvalues E n of 

H N \n) = E n {N)\n) (1.1) 

where |n) stands for an eigenvector of the Hamiltonian H N . 

Now the condition of thermodynamic stability is that for all n 

E n > -AN (1.2) 

where A is a positive constant and indeed, many-body physics is only interesting for 
those systems where the ground state energy per particle 

e = Urn Eqs/N (1.3) 

exists. Thus to consider the spectrum of (1.1) in the N — > oo limit we are most 
interested in 

Urn (E n - E GS ). (1.4) 

In the examples above of the free gases, Landau theory, and BCS theory there 
are an infinite number of levels E n for which the limit (1.4) exists. These are what 
we call the low-lying energy levels. However for these and most other translationally 
invariant systems the much stronger property holds, that 



Urn (E n -E GS ) = V ei(Pi) (1.5a) 

N—>oo A — ' 
i,rules 

and 

Hm(P n ~P GS )= V P t (1.56) 

N^oo z — ' 
i,rules 

where e»(P) is a set of single particle energy levels, P is a momentum, and £\ rules 
indicates that the energy levels are additively composed with a certain set of rules, 
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which embody bosonic or fermionic restrictions and internal quantum numbers such 
as spin or isospin. This additive composition of the order one terms in the energy 
levels in the N — » oo limit and the relation of energy to momentum of (1.5) is what we 
mean by a "quasi-particle" spectrum. (The phrase "quasi-particle" has many related 
usages in many-body physics. In this paper we will use only definition (1.5)) 

The property of possessing a quasi-particle spectrum occurs in a very large number 
of translationally invariant many-body systems. It occurs in free gases, is a key 
ingredient in the Landau theory of the Fermi liquid, and is an important consequence 
of BCS theory. However, regardless of how ubiquitous its occurrence in many-body 
physics the quasi-particle spectrum does not follow from general properties such as 
hermiticity and a finite range assumption on the potential alone. Indeed, in BCS 
theory it can be argued that the quasi-particle nature of the spectrum is almost an 
accident because there seems to be no a priori reason that the energy of an excited 
pair should be precisely the sum of the energies of 2 quasi-particles having opposite 
momentum. The question is thus raised of whether there exist many-body systems 
whose spectra fail to have the quasi-particle form. 

The first example of a non quasi-particle spectrum was discovered in the /V-state 
superintegrable chiral Potts model. This model was originally introduced to study the 
quantum commensurate-incommensurate phase transition [1] and as a representation 
[2] of Onsager's algebra [3]. Its spectrum was studied by recursion relations [4] and 
by functional equations [5,6,7] and this led to the discovery of the non quasi-particle 
excitations. Subsequently one of the most important features of the non quasi-particle 
spectrum was demonstrated to follow from Onsager's algebra alone [8]. However 
the complete spectrum has not been presented in detail because of the lack of a 
completeness relation that counts all states. The purpose of this paper is to present 
this completeness relation for the three state case and to discuss the physics of the 
non quasi-particle excitations. 

We define the model and summarize the formalism needed for the computation 
of the spectrum in section 2. In section 3 we present the rules which count the states, 
and in the appendix demonstrate that these rules are complete. These counting rules, 
indeed, have a very interesting relation to the 5 = 1 XXZ model [9-15] with anisotropy 
parameter 7 = 7r/3 which we discuss in section 4. In section 5 we compute all the 
low-lying excitations in the massive phase. In particular, the complete result which 
includes the non quasi-particle excitations is given in (5.13). However, even though 
the non quasi-particle excitations are a new phenomenon they can be regarded as a 
generalization of the excitation spectrum of the BCS system [16]. This interpretation 
is discussed in section 6. We conclude in section 7 with a summary of the new 
phenomena found in this model. 
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2. Formulation 



The superintegrable chiral Potts chain is the specialization of the general inte- 
grable chiral Potts chain 



M N-l 

H = "EE {«n(^-) n + a n (ZjZ} +i r} (2.1) 

j=l n=l 

with 

= I N <g> • • • <g> X J ' th <g> • • • <g> 7jv (2.2a) 

and 

Zj = I N <g> • • • <g> Z Jt/l ® • • • ® 7jv (2.26) 

where 7^ is the N x N identity matrix, the elements of the N x N matrices Z and X 
are: 

Z/,m = <Si,m^~\ (2.3a) 

- <*i,m+i (mod N), (2.36) 

w - e 2m ' N (2.4) 



and 

a„ = exp [i(2n - N)<f>/N]/ sin {nn/N), (2.5a) 

a n = A exp [i(2n - N)4>/N}/ sin {nn/N), (2.56) 

cos^ = Xcos(f> (2.6) 

to the special point 

^ = ^ = -nj2. (2.7) 

The formalism needed to compute the eigenvalues of (2.1) in the superintegrable 
case (2.7) has been presented extensively in ref. [7]. We summarize here only the 
results needed to study the completeness of the spectrum. For detailed derivations 
the reader is referred to ref. [7]. 

The chiral Potts chain is solvable because of its relation to the 2 dimensional 
chiral Potts statistical mechanical model [17-21] whose transfer matrix is 

M 

where the local Boltzmann weights 

Wp, g ( n ) _ tt fd p b q - a p C q LO j 



= T]( 77 ^ ) (2.9a) 

± ^ \Opd q — C p a q Lj3 ) 



3=1 
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and 

W p q (n) / ujcipdq — d. 



TTf^-VV-n {29b) 
1\\ c p b q -b p c q uji J 



W p h q (0) , 

satisfy the star-triangle equations. Here a p ,b p ,c p ,d p and a q ,b q ,c q ,d q lie on the general- 
ized elliptic curve 

a N + Xb N = X'b N , \a N + b N = \'c N (2.10) 

with 

A' = (l-A 2 ) 1 / 2 . (2.11) 

These Boltzmann weights depend only on the difference lj - I'j (and satisfy 
W v ' h (n+N) = W v ' h (n)) and the Hamiltonian (2.1) is invariant under the global rotation 
of all spins by 2tt/N. Thus the eigenvalues of T and H may be classified by the Z N 
charge eigenvalue e 27riQ/Ar . We further note that (2.1) is obtained from (2.8) in the 
limit q — » p as 



, JV-1 



T « = + 2 " M &" /2 1 + uH + ° ( " 2 » <2 ' 12) 

and that 

e 2id>/N = _l/2 a P c P ^ c 2i4>/N = _l/2 a P d P _ , 2 .13) 

b p d p b p c p 

The method of computation of the spectrum of (2.1) for N = 3 presented in ref. [7] 
is based on the fact that T p , q satisfies a functional equation. This equation is derived 
in detail in ref. [22]. For the superintegrable case (2.7) 

a p = b p , c p = dp (2-14) 

and 



Cp 



(l-A)]l/2iV 



(1 + A) 



(2.15) 



In this case we explicitly factor T p , q into its zeros and poles using 

T - ytd tl t n (2 16) 

p ' q ~ [(f^) 3 - l]M P-i l z - iD J 

with 



^=[(l + A)/(l-A)] 1 /e (2.17) 

(where the index q on a, b, c, d is suppressed) and the functional equation is explicitly 
written as 

T p N g Tp N Rq T p % q = * M e- iP {{f/ 1) M (^ 2 - 2 - V" 7 ?* 

+ (^ V W 1) M (^ 2 - - l) M T p N R2q (2.18) 

+ ( ^ 2 - 1)M( ^ 2w - 1)MT ^ } ' 
5 



where R is the automorphism 

R(a, b, c, d) = (b, uja, d, c) (2.19) 
and P is the momentum which is derived from T p , q as 

e~ lP = limTp Rq . (2.20) 

The functional equation (2.18) is solved by the ansatz (note that the sign of vi 
here is defined opposite to that of ref. [7]) 

3 m P -i _ 2 ab rn E 



,».3« ( ^ ( ^ ( ^n ( i^)ri(^)-{^ ± ^} ,, 21 , 



where the ± signs are chosen independently for each of the terms in the product, and 
P a ,Pb,P c are positive integers. The numbers v t satisfy 



(^—^) M = - w -(^) TTf "*-^ ) . (2.22) 
The wi are obtained from 



\;\ - A)" + — 

</ 

where the U are the roots of the polynomial 



wf = - A (l-\y + T - ¥ (2.23) 



P(t) = t-^+^i^H - l) M (ut - l) M w p « +p > f[ Q_ f 3 ) 



+(*-i) M (^-i) M n(rr^i) ^ 

m P 1 2+ 

+ (t-i)«( W t-i)" W -^)n(^^)}. 

The eigenvalues of H are thus finally obtained as 

m E 

E = A + BX + 6^2± Wl (2.25) 
1=1 

where 

A = 3(2P c + m E ) - 2M (2.26a) 
B = 2(P b — P a ) - A (2.266) 
6 



and the ± signs are all independently chosen. In addition the momentum P is obtained 

as 

-l - cj 2 vr 

LUVl 



The ansatz (2.21) gives more possible eigenvalues than the allowed number of 3 M . 
The completeness problem is the determination of those values of P a , Pb, P c , m p , m E and 
those possible solutions vi of (2.22) which actually give eigenvalues of H and T. 

It should also be noted that if in (2.22) we set 



v k = e Afc (2.28) 



we obtain 



\sinh\(\ k -lf)) " ll sm ^i( Afe -A ; -^f)" 

Up to the phase factor in front this equation is the Bethe's equation for the spin 5 = 1, 
anisotropy 7 = n/3 XXZ chain with M sites and \S Z \ = M - m p [12-15]. In the case 
where M-m p -P a -P b = (mod 3) eqn. (2.29) is identical with the 5 = 1 equation with 
periodic boundary conditions Thus the study of completeness of the superintegrable 
chiral Potts model is closely related to the completeness of this XXZ chain. 
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3. Completeness Rules 

It is a feature common to all studies of eigenvalues which use functional equations 
that further information beyond the functional equation is needed to determine the 
multiplicities of the eigenvalues (and in particular to determine whether or not the 
multiplicity is zero). It is to be expected that such additional information can be found 
by studying the eigenvectors accompanying the eigenvalues. However the production 
of the eigenvectors, though of great importance, is also a problem of substantially 
increased complexity and in practice it is desirable to obtain the completeness rules 
in advance of an eigenvector computation. 

Such a study of completeness rules was recently done for the non-chiral limit 
= = of (2.1) [23] by combining a finite size study of eigenvalues of chains of 
size up to M=7 with a generalization of the string hypothesis [24-27] approach to 
the counting of solutions of Bethe's ansatz equations. We will here follow a similar 
approach for the superintegrable case (2.7). First we will use the results of a finite 
size study to obtain a set of rules for the allowed integers P a , P b , P c , m E and m p , then we 
will use a finite size study of the v t to find a set of rules which counts all eigenvalues. 
In the appendix we show that these rules derived from systems of size M = 3,4,5,6 
sum up to the total number of states 3 M for all M. 



A. Rules for P a , Pb, P c , hie, and m p 

From (2.26a) we see that instead of specifying P c it is sufficient to specify the 
coefficient A in (2.25). We thus summarize the results of the finite size study of ref. 
[7] by presenting in table 1 the allowed values of A + BX,P a , and P b . For each set the 
integers, m E are either even or odd. This parity is also listed in table 1. We note that 
within each set there is a sum rule on the quantity 

3m E + 2m p + P a + P b . (3.1) 

The value of this sum rule is also listed in table 1. In addition we include the value 
of M - P a - P b - m p which determines the phase factor in (2.29). 



B. Rules for vi 

In table 2 we give for M = 3, 4, 5 and 6 the number of allowed sets of v k (denoted 
by N) which correspond to the allowed values of P a ,P b ,m E , and m p given in table 1. 
From (2.21) each set of v k corresponds to 2 mE eigenvalues of T so that the N sets 
give 2 mE x N eigenvalues. The numbers N are found from table 1 of ref. [7]. These 
numbers N must count the number of allowed solutions of (2.22). 
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To study the solutions of (2.22) we have computed the zeroes of T by means of 
the technique of ref. [23] for M = 3,4, 5 and 6. From that study we find that even for 
M=3 there are only 3 classes of solutions: 

1. real positive solutions (denoted by v£) 

v k >0 (3.2a) 

2. real negative solutions (denoted by v^) 

v k < (3.26) 

3. complex conjugate pair solutions (denoted by v\\ s ) 

Argv k ~ ±7r/3. (3.2c) 

Following the conventions of ref. [23] we refer to these as positive parity, negative 
parity, and (positive parity) 2 strings, and denote the number of such roots in a given 
eigenvalue as to + ,to_, and m 2s respectively. Clearly 



TO + + TO_ + 2to,2 s 



(3.3) 



Following the practice of refs. [12-15] and [23-27] we make these three classes of 



roots explicit by rewriting (2.22) for v k and v k as 



(4^r=(-D 

V V, — LJ / 



1 \ M 



Pa-Pb-rrip , 



n 

i=i 

rri2s 



iii_ _i_ 

n( v kU~Vi \ 
Avrw-vf) 



i -i 
± 



k 
2s* . 



m 2s + 2s 9s* 



Also, multiplying the equations for vl s and w| s * together we obtain 

- 1- 



2s* 



2s 



M 



2(M-P a -P b -m p ) . 



rp 

m,2s 

n[( 



,2s 



,2s*, 



2s — 2s* 



w — w. 



,,2s* 



vl s UJ 



vI s )\v? s lj - vl s *J \v? s *uj 



,,2s* 



,2s 



,2s* 



.,2s* 



(3.4a) 



(3.46) 



We continue to follow [12-15] and [23-27] by taking the logarithm of (3.4). This 
introduces integers or half integers [28] as the various branches of In 1 or ln(-l). These 
counting integers are of great importance in classifying the solutions. Their properties 



are dependent on the choices of branches for logarithms. The principle used to choose 
branches is discussed in reference [23]. Thus we define all logarithms to obey 



and we use the definitions 



-77 < Im In x < 7r (and In 1 = 0) 



, t+(t , + )= / „(-_ r _) > 
i>~u; — 1 



zi_(v ) = ln(— 

/ 2jn , IV^fc^ _ 1 V w fc S * w _ ^ 



and 



iQ +,+( v k' v i) = ln 



i@-,-( v k ' v , 
iQ+,2s{vt^ s ) = ln 



i&-,2s{v k ,v? s ) = In 



In 



V.UJ- vt 



Vu U) — V 



2,s* 

I 



i®2s,2s{v 2 k s ,v^ s ) = In 



2s* 



( v k u ~ v i \ / v k u-v,' V 
Vu 2s ui - v k ' \vf s *u> — ) - 

s vl s uj-v? s \ s v 2 k s *uj-v? s \ 

V W ; 2S W — U| S / V W; 2S W — V k S * J 

/ v£u--v^_ \ / vl s *oj-v? s * \ 
\v? s *ll> — v1 s J \v? s *lu — v? s * J 



and 



(3.5) 

(3.6a) 
(3.66) 
(3.6c) 

(3.7a) 
(3.76) 
(3.7c) 
(3.7d) 
(3.7e) 

(3-7/) 
(3.8) 



8 m ,„(j;,!/) = -9„,m(j/,ar)- 
We note in particular that the minus signs in these definitions are chosen so that 

= = t 2s (-^ ±1 ) = 0. (3.9) 

Taking the logarithms of (3.4), we define the counting integers P k by 

Mt+(v+) = i lnLU M - p «- Pb - m ? + 2tt/+ 

m + m- m 23 

+ E +.+K + '^ + ) + E e +.-K + 'V) + E0+,2sK + ,^), 



(3.10a) 



1=1 



1=1 
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Mt_(«fc ) = j lnu M - p «- p »- m r + 27r/ fe -+ 



E -,+K~>O + E e -,-K"^D + E e-, 2a («fc ,«?*), 
1=1 1=1 1=1 

m + 

Mt 2s (v 2 k s ) = - inJ*W-Pa-i%-m,) + 27r / fe 2s + J2 © 2s ,+(^ s , 
2 l=i 

m- m-2s 

+ E e 2s ,-(^>r) + E ^, 2s (^, «?•), 



(3.106) 



(3.10c) 



i=i 



where 



T+ . ran integer "i . t „, , .. , • reven^i ,„,, . 

J * 18 la half-integer} ^ ^ + 1 + - + ^ | odd |. (3.11a) 

.._ . ran integer l - f , , • f even "i ,„ , , , . 

J * 18 la half-integer } lf m ~ + m >° + 1 18 { odd }" ( 3 ' 115 ) 

r9s . r an integer 1 , • f even ~i /„,,•> 

7 " 18 la half-integer} ^ ™- + ™». + 1 is { odd }. (3.11c) 

We are now able to discuss the rules for completeness of the vi. As will become 
apparent, the cases Q = 0, 1 and 2 must be treated separately. 



1. Q=0 

The completeness rules are simplest for Q=0. Following the procedure of ref. 
[23-25], we count solutions of (3.10) by finding the allowed values P k and assuming 
a monotonic relation between P k and v\. We thus find the limits within which the 
integers P k can lie by examining the extreme cases of (3.10), where \v\\ = and oo. 
This is easiest for M = m p = (mod 3), P a = P b = 0, where we find that the number of 
allowed values for I k ,I k ,I k s are 

iV + (m + ,m_,m 2s )- 3V pJ , (3.12a) 

V m + / 

AT t \ fh( 2M + m p) - m 2s\ ,,. n 

N-(m + ,m-,m 2s ) = [ 6 F , (3.126) 



where 



,\{2M + m p )-m + -m 2s , 
iV 2s (m + ,m_,m 2s ) = ( J ^ ), (3.12c) 



(3.13) 



y bJ (a -&)!&! 

Thus the total number of solutions to (3.10) with J^, I k , I k s separately distinct 

is 

7V(m + ,TO_,TO 2s ) = N+N-N 2s . (3.14) 
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This counting, however, does not give the correct result. This can be seen for the 
cases M = 6, nip = 3 and 6 where from (3.14) we find the total number of states 
satisfying the sum rule (3.3) is 50 and 141, whereas we see in table 2 that the correct 
number of states is 46 and 43. Thus some further correlation among the integers P k 
is required. Guided by ref. [23], we make the additional conjecture 



If s + -II 



(3.15) 



We compute I 3 k from (3.10) for the solutions v k obtained from our M = 6 finite size 
study and find that indeed (3.15) is satisfied. Furthermore, we find that the counting 
of states given by (3.12) and (3.15) does indeed give 46 for m p = 3 and 43 for m p = 6 
as required. 

For the case M = 1,2 (mod 3) the rules are almost as simple, except we must use 
the function [x] which is the greatest integer less than or equal to x. Thus we obtain 
the result valid for Q = and all M, that for all states allowed by table 1, the number 
of allowed solutions to the Bethe's equations (3.4) for v t are 



N(m + ,m-,m 2s ) = 



M+l-m p 



2M+m„ 



- m + - m 2l 



2M+m v 



(3.16) 



m+ — m_ — TO2., 
m 2s 



where I k and /| s obey (3.15). 



2. Q=l 



For Q=l our procedure is very similar, except now, as may be inferred from table 
1, several separate cases must be considered. 



a. M = 2 (mod 3), P„ = P h = for all m p 

Here the counting rule is identical with (3.16): 



N 



M+l-TOp 

3 

m + 

~ 2M+m p 
3 



2M- 



771+ - 7712s 



777- 

777 + — 777_ — 7712; 

m 2s 



(3.17) 
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b. M = (mod 3), P„ = 2,P b = 0,w p = 0,1 (mod 3) or 
M = 1 (mod 3), P„. = 0, P h = 2, m r = 1, 2 (mod 3) 



N(m + ,m-,m 2s ) 



M-m 



+ 171- 
2M-2+m p 



2M~2+m p 
3 



- m + - m 2l 



— m + — m_ — m 2s 
m 2s 



(3.18a) 



a M = (mod 3), P„. = 0^= l,m p = 0,2 (mod 3) or 
M = 1 (mod 3). P„. = l.Ph = Q.m v = 0.1 (mod 3) 



N(m+, m-,m 2s ) 



+ m_ 



2M-2+m p 
3 



2M-2+H1. 



- m + - m, 2 t 



- m + — m_ 
m 2s 



™2 ; 



(3.186) 



3. Q=2 



The rules for Q=2 are obtained by an analogous procedure. The results are: 



a. M =. 1 (mod 3), P„ = P h = for all m p 



N(m + ,m- 1 m 2s ) 



M— m p +l 
3 

2M+n 



2MH 



m + — m 2i . 



m + — to_ — m 2s 
m 2s 



(3.19) 
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b. M = (mod 3), P„ =0,P h = 2,m p = 1,0 (mod 3) or 
M = 2(mod 3), P n = 2,P h = 0,m r = 0,2 (mod 3) 



N(m + ,m-,m 2s ) 



M-m 



+ 171- 
2M+m„-2 



2M+m„-2 



- m + - m 2; 



— m + — m_ — m 2s 
m 2s 



(3.20) 



c. M = (mod 3), P„ = 1^= 0,m p = 0,2 (mod 3) or 
M ee 2 (mod 3), P, = 0, P, = 1, m p = 1, 2 (mod3) 



N(m+, m_,m 2s ) 



M-m v +2 



+ m_ 



2M+m p -2 
3 



2M+ro„-2 



- m + - m 2 < 



- m + — m_ — m 2; 
m 2s 



In all cases the exclusion rule (3.15) holds. 

In the appendix we show for arbitrary M and Q that 

2 mE N(m+,m-,m 2s )=3 M - 1 . 

m+ ,m- ,ni2 s 

Thus these rules correctly count the number of states for all M. 



(3.21) 



(3.22) 
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4. Completeness for the S =1, "/ = 7r/3, XXZ spin chain 



In section 2 we remarked that the Bethe's equation (2.29) of the 3 state superinte- 
grable chiral Potts model is (up to a possible phase factor) identical with the Bethe's 
equation which arises in the solution of the 5=1 anisotropic XXZ chain with 7 = n/3. 
Thus it is to be expected that there is a close connection between these two models. 

In the superintegrable chiral Potts model the eigenvalues are grouped into sets of 
size 2 mE . Each eigenvalue of the set is specified by the same solution of the Bethe's 
equation (2.29) and they differ by the 2 mE independent choices of the m E ± signs in 
(2.21) and (2.25). We will see in this section that a completely analagous phenomenon 
occurs in the XXZ chain with 7 = n/3. Namely, to each solution of the Bethe's 
equation (2.29) there also corresponds a set of eigenvalues for the XXZ hamiltonian. 
However, in distinction with the chiral Potts chain, the eigenvalues in these sets for 
the XXZ chain are degenerate. 

The 5=1 anisotropic integrable spin chain is defined by the hamiltonian [9] 



M 

H = ~y^{Sj ■ S j+1 - (Sj ■ S j+ iY 



3=1 

-tz\2 



-2 8irSi[S?Sj +1 - (S*S* +1 Y + 2(Stf] (4.1) 
2( COS7 - 1)[(5|5| +1 + S]S] +1 )S*S* +1 



, QZ QZ ( QX QX , QV QV \] \ 



where periodic boundary conditions are imposed and the single site 3x3 matrices 

gx, y ,z are 



(4.2a) 

S y = — I 1 1 I , (4.26) 




and 

S z = I I . (4.2c) 



The operator 



M 



s z - J2 S 3 ( 4 - 3 ) 



commutes with (4.1) because of rotational invariance in the XY plane. The eigen- 
values of this spin chain have been extensively studied [12-15] and in particular it is 
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shown in refs. [12] and [13] that the energy eigenvalues are given by 



M-\S'\ 

E = sin 2 (2 1 ) V =- =- (4.4) 

~i sinh\{Xj + 2i'))sinh\{\ j - 2i-y) 

and the corresponding momentum eigenvalue is given by 

sinh^X j -2iy) 
where the X k satisfy the Bethe's equation 

sinh^{X k + 2vyh m = 1 s inh \{X k - A, + 2z 7 ) 

sinh\{X k -2i 1 )\ ~ J! sinh ^x k -X l -2i 1 )' 

When 7 = 7r/3 this Bethe's equation is identical with the equation (2.29) of the super- 
integrable chiral Potts model for the cases where M - P a - P b - m p = (mod 3). 

The spectrum and thermodynamics of (4.1) have been extensively studied in refs. 
[12-15] and for 7 = the question of completeness has been solved in ref. [29, 30]. 
When j/tt is rational a set of degeneracies between energy levels occurs which can 
be studied by means of the quantum algebra U q {SU(2)\ at a root of unity [31,32]. 
However, for our purposes here of comparing the XXZ chain with the superintegrable 
chiral Potts model we will not use this general formalism and will content ourselves 
with the examination of a six site chain. 

We consider first those cases where the phase factor in equation (2.29) is unity. 
From table 1 we see that this occurs in Q = for m p = 0, 3,and 6; in Q = 1 for P a = 2, 
P b = 0, m p = 1, 4 and for P a = 0, P b = 1, m p = 2, 5 and; in Q = 2 for P a = 0, P b = 2, m p = 1, 4 
and for P a = 1, P b = 0, m p = 2, 5. Thus all values of m p = M - \S Z \ occur. 

To compare the number of states in the XXZ chain for given S z and P with the 
number of allowed solutions of the corresponding superintegrable chiral Potts chain 
we must relate m p to S z and P C p to P s =i- The correspondence m p = M - \S Z \ is valid 
for all m p . To find the correspondence of momentum we compare the momentum 
formula of chiral Potts (2.27) with the momentum formula of the S=l chain (4.5) to 
obtain 

e iPcp = e iPs=1 LO m "~ Pb . (4.7) 
Thus we find for all cases where M - P a - P b - m p = (mod 3) that 

if to p ee (mod 3) then P s =i = P C p (4.8a) 

and 

2ir 

if to p = 1, 2 (mod 3) then P s =i = P C p - -rr- (4.86) 
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Using these correspondences we list in table 3 the number of states of the 5 = 1 model 
and the number of allowed solutions of (2.29) for the superintegrable chiral Potts 
model obtained from table 1 of ref. [7]. In terms of the Ps=i defined in (4.8) the 
number of chiral Potts states is symmetric in ±Ps=i even though for m p = 1, 2 (mod3) 
they are not symmetric in P C p. 

It is obvious from table 3 that the number of eigenvalues of the 5 = 1,7 = n/3 
XXZ chain is significantly larger that the number of solutions of the Bethe's equation 
(2.29) allowed by the counting rules of sec. 3. The question is thus to find where the 
extra solutions of the XXZ chain come from. 

In table 4 we compute all eigenvalues of the XXZ chain (4.1) with 7 = ir/3 for 
a chain of six sites. We have also computed that subset of eigenvalues obtained by 
substituting in the energy formula (4.4) the solutions of (2.29) with phase factor one 
which exist for the superintegrable chiral Potts chain. This subset of eigenvalues is 
underlined in table 4. 

For further discussion it is necessary to consider \S Z \ =0,3 and 6 separately from 
\S Z \ = 1,2,4, and 5. 

First consider \S Z \ = 0,3, and 6. Inspection of table 4 shows that the single under- 
lined eigenvalue obtained in P = 0, S z = 6 from the superintegrable chiral Potts model 
recurs 4 times in S z = 3, 6 times in S z = 0, 3 times in S z = -3 and once in S z = -6 for 
a total multiplicity of 2 4 . Similarly the underlined eigenvalues which occur in S z = 3 
recur twice in S z = and once in S z = -3 for a total multiplicity of 2 2 . This is an 
illustration of the rule that an eigenvalue obtained from the chiral Potts chain which 
occurs singly in the sector S^ ax will recur in all sectors S z congruent to mod 3 
with a total multiplicity 

This is precisely the factor 2 mE which occurs in the superintegrable completeness sum 
(3.22) and hence the computation of section 3 shows that the number of levels in the 
sector S z = 0(mod3) is 3 M -\ 

From the point of view of the Bethe's equation (4.6) these 2 2S ™-/ 3 degenerate 
eigenvalues correspond to solutions obtained by letting some of the go to ±00. 
This mechanism of multiple occupancy of A fe = ±00 is the same phenomenon that 
occurs in the isotropic case 7 = for the cases 5=1 [29,30] and 5 = 1/2 [24,26,27] and 
in the Hubbard model [33]. 

For the case S z = l,2(mod3) the situation is slightly more complicated. Consider 
first S z = 5. Here table 4 tells us that in each of P = 0,±27r/6 and tt there is a single 
level which comes from the superintegrable chiral Potts model and that this level 
recurs twice in S z = 2 and once in S z = -1 (and similarly the states in S z = -5 recur 
twice in S z = -2 and once in S z = l). Thus the total multiplicity is 2 2 . Thus we have 
an identification of all the chiral Potts states in Q = 1 and 2 and m p = l,2(mod3) with 
the states in 5 = 1. 

We have now exhausted all the chiral Potts states with phase factor one. However 
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not all the states of the S = 1 model have been been accounted for. To obtain 
the remaining levels we must consider the remaining solutions of the chiral Potts 
equation (2.29) where the phase factor is w ±:L (see table 1). We see in table 4 that 
the eigenvalue 18 in P = ±2ir/3 for S z = 5(-5) recurs 3 times in S z = 2(-2), 3 times 
in S z = -1(1) and once in S z = -4(4) for a total multiplicity of 2 3 . This is identical 
with the multiplicity of the eigenvalues of table 2 with Q = 1, P a = 2, = 0, m p = and 
Q = 2, P a = 0, P b = 2, m p = 0. The remaining eigenvalues are degenerate in S z = 2(-l) 
and -1(1) and are obtained by substituting the chiral Potts values of A fc into the 
energy formula (4.4) . Thus we have demonstrated for M — 6 that all the eigenvalues 
of the 5 = 1 chain with 7 = tt/3 can be obtained from the same solutions of the 
Bethe's equation (2.29) that gave the superintegrable chiral Potts eigenvalues. The 
degeneracies in the 5 = 1 chain correspond to multiple occupation of A fe = ±00 and 
the same exclusion rule (3.15) found for chiral Potts holds also in the XX Z chain. 
This phenomena holds for all M. In general the 2 mE states which correspond to the 
choice of the m E ± signs in the chiral Potts model give in the 5 = 1 model degenerate 
eigenvalues for values of S z which are congruent mod 3 and in each S z sector the 
degeneracy is given by a binomial coefficient. 
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5. Excitations in the massive phase 

In refs. [6,7] it was shown that the superintegrable chiral Potts chain (2.1) has a 
massive phase if 

< A < A c = .9013... and 1 < A < A" 1 . (5.1) 

In this phase, the spectrum of excitations for the sector Q = is computed in detail. 
In this section, we combine the computations of ref. [7] with the completeness rules 
of section 3 to explore the physics of the model. 

We begin by recalling that if for Q = we consider the excitations in (2.25) where 
all minus signs are chosen, then the result (4.9) of ref. [7] for the energy eigenvalues 
is 



Urn {E(P,X,m p ) - E GS } = 2m p \l - A| 

M— >oc 



„ „|i±A| 2 /3 m p „ 



x 



(5.2) 



4 \-d-A)f /2 



and the corresponding momentum is 

«-"-n(^) M 



(where we restrict our attention to M = (mod 3) so that m p = (mod 3) and 
P a = P b = 0). 

To make contact with the completeness discussion of section 3, we consider v+ 
and v~ together as 

v = v r where - oo < v r < oo (5.4a) 

and the 2-strings as 

v = v 2s e ±7ri/3 where v 2s > (5.46) 



and define 



| l + A ,2/3 

3 f ll - xl , f uv r lo z v t 1 

s r = 2 l - A + - / dt\ - + x 

111 - A| + - / „, , 



and 



(5.5a) 



3 I-It^I 27 ' 5 W 2-[4( v 2. t) 2_ w 2. t+1 ] 

e 2s =4 1-A +- / dt- 1 V ' 



4A 



t 3 - 1 



(1-A) 



1/2 



(5.56) 
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Thus we rewrite (5.2) as 

m r mis 

Jim {E(P, A, m r , m 2s ) - E GS } = Y] e r (w, r ) + V] e 2s (vf s ) (5.6) 

1=1 i=i 

where 

m r = m + + to_ (5.7) 

and the v\ and vf" satisfy the exclusion rules vf ^ v r k and vf s ^ vl s . Similarly, we 
rewrite (5.3) by defining 



and 



_ ip 2s _ / 1 - w 2 eW 3 w 2s x / 1 _ u 2 e -^/3 v 2s 
_ V 1 _ we"/^ 28 / VI - w e -"/ 3 i; 2s J 



(5.86) 



Thus (5.3) becomes 



p = E p < r + E p ' S (mod 2^) (5.9) 



and we note from (5.4) and (5.8) that 

< P r < 2tt (5.10a) 

whereas 

y < P 2s < 2tt. (5.106) 

When (5.1) holds, both e r and e 2s are positive. The phase transition at Xf 1 occurs 
because the gap in e r (but not in e 2s ) vanishes [6,7]. 

To complete the presentation of the energy spectrum of (2.1) it remains to consider 
replacing an arbitrary number of minus signs in (2.25) by plus signs. This is equivalent 
to adding 2wi to the energy where wi is given by (2.23) and t t are the roots of the 
polynomial (2.24) in the ground state m p = 0. In ref. [7] it is shown that in the 
ground state asM^oo the roots are such that tf is less than minus one. Thus we set 
-tf = (vf) 3 and define 



4A 



1/2 



eM = 6 [iTW + {1 ~ X) ^ ' (5 ' U) 

with 

1 < vf < oo (5.12) 

where we note that e c is positive for A ^ 1. Then, denoting the number of changed 

minus signs as m c , we obtain the complete explicit formula for the energy and mo- 
mentum of all excitations as 

Urn {E ex (P) - E GS } = V e r (vl) + V e 2s (vf s ) + V e c (vf) (5.13a) 
1=1 1=1 1=1 
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and 



Urn P = Y, P l + T, P i 

1=1 1=1 



,2s 



(5.136) 



where 



— oo < v r < oo, < v 



,2s 



1 < « c , 



(5.13c) 



m r + 2m2 S = 0(morf3), 



(5.13d) 



e r , e 2s , e c , P r and P 2s are given by (5.5), (5.11) and (5.8), and the vi satisfy the 
exclusion rules vf ^ t>£, v 2s 7^ w 2s and ^ v%. In addition, we note that corresponding 
to the exclusion rule 1^ ^ -If s (3.15), there will be an exclusion rule between negative 
v r and v 2s of 



We are now in a position to discuss the physics of the excitation spectrum of 
the superintegrable chiral Potts model. In particular, we return to the discussion of 
the introduction and compare the excitation spectrum (5.13) with the quasi-particle 
spectrum (1.5). 

Consider first the special case m c = 0. Then (5.13) has precisely the form (1.5) 
where the excitation energies e r (P) and e 2s (P) are obtained by eliminating v r and v 2s 
between (5.5) and (5.8). Thus, using the definition of "quasi-particle" given in the 
introduction, we may say that e r (P) and e 2s (P) are the energies of single quasi-particle 
excitations . 

However, there are two striking differences between these quasi-particles and all 
the previously discussed quasi-particle excitations; namely, the momentum of the 2s 
excitations is only allowed from (5.10b) to be in the range 27r/3 < P 2s < 2ir and from 
(5.14) there is an exclusion rule between real excitations of momentum < P r < An/3 
and 2s excitations. It thus seems completely fair to say that even though real and 2 
string excitations have the quasi-particle form (1.5) for energies and momentum, the 
rules of combination are qualitatively different from those derived from the intuition 
that is typically implied by the word "quasi-particles." 

But even more striking than these new combination rules for energy levels is the 
form that the energy levels take when m c ^ 0. In this case, we have a contribution e c (v c ) 
to the energy which depends on one parameter but makes no contribution whatsoever 
to the momentum. This violates the quasi-particle form (1.5) and is what we mean 
by non quasi-particle excitations. 




(5.14) 
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6. Non Quasi-Particle excitations and BCS Theory 

The superintegrable chiral Potts model is the first system in which non quasi- 
particle excitations, as defined in section 5, have been observed in either exact or ap- 
proximate computations. Nevertheless these excitations have a striking resemblance 
to the Cooper pair excitations of the BCS theory of superconductivity [16]. 

In section 3 of their famous 1957 paper [16] Bardeen, Cooper, and Schrieffer 
classify the states of a superconducting system into three types: single particle, ground 
pairs, and excited pairs. The implication of this classification scheme is that single 
particle states and Cooper pairs are logically distinct. 

Single particle states are specified by a momentum k. Excited Cooper pairs are 
specified by a vector k'. However the momentum of the Cooper pair is zero. It is a 
most important result of computation that in equation (2.52) of ref. [16] "the energy 
to form an excited pair in state fc"' is 2Ep where E^, is the energy of a single particle 
excitation given by equation (2.50) of ref [16]. The result of this computation is 
mandatory for the interpretation of a Cooper pair as two quasi-particles with opposite 
momenta. This is the interpretation which allows the BCS spectrum to be interpreted 
in the quasi-particle form (1.5). 

But the clear reading of ref. [16] is that this identification of the energy of 
the excited Cooper pair with 2E%,is a result of computation, and does not follow 
from model independent general principles. The similar situation would arise in the 
superintegrable chiral Potts model if 

2e r (P) = e c (P). (6.1) 

We have seen in section 5 that this relation does not hold for the 3 state super- 
integrable chiral Potts model. However, if we consider instead the two state case 
where the chiral Potts model reduces to the Ising model the relation (6.1) does in- 
deed hold and, as is well known [3, 34, 35], the Ising model spectrum does indeed 
have the quasi-particle form (1.5). Indeed, it is significant that in ref. [35] Schulz, 
Mattis and Lieb note the strong similarity of the Ising model with the classic works 
on superconductivity [16, 36-38]. 

The results of this paper, of course, stand by themselves without reference to 
BCS theory. Nevertheless the above discussion of the presentation of reference [16] 
suggests that the excitations e c can be thought of as the generalization of Cooper 
pairs to a situation where the equality (6.1) fails. 
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7. Conclusion 



In this paper we have discussed three novel properties of the excitation spectrum 
of the superintegrable chiral Potts model: 1) single particle energies (5.5b) which do 
not exist for all ranges of momentum (5.10b), 2) a fermion like exclusion rule (5.14) 
which operates between different single quasi-particle levels and 3) non quasi-particle 
excitations (5.11) which carry energy but make no contribution to the momentum 
(5.13) and cannot be rewritten as the sum of single particle levels. Related phenomena 
certainly exist in the general chiral Potts model (2.1) for 0^0. In addition the 
phenomenon of non quasi-particle excitations continues to exist in the massless phase 
of the chiral Potts model. Indeed, the phase transition at the point A = 1 [7] occurs 
because the gap in the non quasi-particle energy e c vanishes. Thus the chiral Potts 
model can be thought of a providing a new universality class with novel physical effects 
which have been seen first in theoretical studies before having been experimentally 
observed. 

On the other hand there may be reasons which forbid these effects from ever 
appearing in real three dimensional systems. If this is the case we hope that the 
present work will inspire an elucidation of such a non-existence theorem. 
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Appendix 

We shall evaluate the sums (3.22), for all the possible sectors detailed in section 
3. Our method follows ref.[23, 24]. 

Q = 

Here the factor N(m + ,m_,m 2s ) of (3.22) is given by (3.10). 
• M = 0(mod 3),P a = 0,P b = 0,m P = 0(mod 3) . 

Set M = 31, m P = 3p. The corresponding value of m E is obtained from table 1 to be 
2(1 -p). Also, 

r A/[ -I— rn i-i l r M 4-1 — m n l 

l-p, (A.l) 



■2M + mp- 


= 2l+p, 


-M + 1 - m p - 


= I 


l 3 




l 3 





so that the sum (3.22) to be performed for this sector is written as 

'21 +p — m + — m 2s 



« (o = E E 22( " p f 

«=fl m. m. i to.i_ 



TO_ 



X 



p—0 m_ ,m_)_,m2 

'21 + p — m + — m 2s — to_^ fl—p + m 
m 2 



to+ 



= E 22( " P) E 

p—0 l v m_,m+,m2s 



5(2m 2s + to_ + to + - 3p), 

(A2a) 

— 2p + TO_ + TO 2s \ /2Z — 2p + TO 2s 



TO 2;S 



^ — 2? — I - 777, 

<\ ^ ~ ) (5(2m 2;s + to- + to + - 3p) }>, 



(A26) 



where 5(a) = 5 O)0 . We use the relation 



21 - 2p + m_ + m 2s \ /2Z - 2p + m 2s \ _ (21 - 2p + m_ + to 2s \ (21 - 2p + m. 
m_ / I to 2s J \ m 2s J \ in- 



to rewrite (A. 2) as 

s °(o = E 22( '" p) E ( 

p—0 rr> 



2l-2p + ra- 
ni- 



(A3) 



(A4) 



(2l-2p + m-+m 2s \ fl-p + m- , , 
> o(2m 2s + m_ + m + - 3p). 

^ \ m 2s ) \ to + 



Then the following two identities, 
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(1 — x 2 ) 2l ~ 2 P +m - +1 



E 

m 2B =0 



2l-2p + m_ + m 2s ^ 
m 2s ' 



and 



l— p+rri- 

(i + x) i - p+m - = 

m_|-— 



I — p + m_ 

m + 



(A5a) 



(A56) 



may be multiplied together, and using the Kroenecker S we get 



p=0 



da; 



2.3P+1 a ,(i-p)+i _ X )2(;- P )+ 



2/ - 2p + to 

to_ 



where the appropriate power of x is picked out by the residue theorem, with a contour 

which surrounds x = and no other singularities. 

Using 



(i - x y 



E 

m_=0 



A - 1 + ra_ 

TO_ 



(1 -2a;) A ' 



(A7) 



(A. 6) reduces to 



da; 



p=0 



= ^ — : </> da; 



^p+i ( X + (i _ 2a: )2(i-p)+i 



2™/ a;(l +a;) i+1 (l - 2a;) 2i + 



p=o 



(1 +x)(l - 2a;) 2 \P 
4^ ) ' 



and the geometric sum is executed to give 



S° (l) = 



da; 
2iri T x 



(1 - 2x) 
(1 - 3x)x 3 



4< +1 a; 3 



(l + a;) i + 1 (l-2x) 2 '+ 1 (l-3x) 



(A8) 



(A9) 



The second term in the square brackets does not have a single pole at x = 0, and 
hence does not contribute to the sum. The first term then gives the desired answer 



• M= l(mod 3),P a = 2,P b = 0,m f = 0(mod 3)orP a = 0,P, = l,m P = 2(mod 3) . 
Here we set M = 31 + 1. For m P = = 3p, P a = 2, P h = 0, the number of states is 



(A10) 



Slm P ^(l) -EE 2 ' 
p=0 m_ ,m+,m2 s 



21 +p — m + — m 2s — m_\ (l — p + m 



m 2s 



(Alia) 
S(2m 2s + to_ + to + — 3p). 
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There are no states for m P = 1 and for m P = 2, P a = 0, P b = 1, the number of states is 



p=lm-.m+.m2 s ^ ' 



X 



p—i m_ ,m+ 

/ 2/ + p — m + — to 2s — to \ (I — p + to_ + 1 
to 2s /V m + 

so that the total number of states is given by 

Si(l) = Si mp =o + Sl, mp =2- 



<5(2m 2s + to_ + m + - 3p + 1), 

(Allc) 



Notice that (A. 11a) is the same as (A. 2). The sums in Si mp = 2 can be performed in 
the same way as S%(1), using relations similar to (A. 3), (A. 5), (A. 7), and the geometric 
sum in (A. 8), to give 

(A12) 
(A13) 



2 • 3 



3Z — 1 



Therefore, 



D l,m P = 2 

S?(0 = 3 3 '- 1 + 2 • 3 3 '- 1 = 3 3i = 3 M ~\ 



• M = 2(mod 3),P a = l,P b = 0,m P = 0(mod 3) or P a = 0,P b = 2,m P = l(mod 3) . 
Set M = 31 + 1. For m P = = 3p, P a = 1, P b = 0, 



2M + m P 



= 2/+p+l, 



M + 1 - TOp 



= / - p+ 1, m E = 2(1 -p) + 1 



while for m P = 1 = 3p + 1, P a = 0, P 6 = 2, 



-2M + top 



= 2l+p+l, 



M + 1 — TOp] 



= 1 - p, top = 2(7 -p). 



(A14a) 



(A146) 



The total number of states is, therefore 

*2S(D = e{ E 2^w( 2 '- p+1 

p=0 \m- ,m + ,m 2 s 



m 2s -m + \ fl - p +1 + m- 

ill J \ TO + 



(21 +p+ 1 - to 2s - to+ - m\ 

o(2m 2s + to_ + to + - 3p) 

V TO 2 S / 



+ 22(i _ p) /2/ +P+ 1 - TO+ - TO 2s \ /2? + P+ 1 - TO 2s - TO+ - TO_ ^ 

\ TO_ J \ m 2s 

I — p + TO. 



5(2to 2s +to_ +to+ - 3p- 1) >. (A16) 

The sum is executed along the same lines as the steps (A. 2) to (A. 10) performed to 
evaluate S$(l). The answer is 



S§(i) = 3 3i+1 = 3 M '\ 



(A.17) 
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Q = i 



• M ee O(mod 3) 

Here, N(m+, m_, m 2s ) is given by (3.18). We set M — 31 and m P = 3p,3p + 1, 3p + 2 
for mp ee 0, 1 and 2 respectively. To evaluate the sum (3.22) for this case there are 
three different terms to be combined. 

(21 +p- l-m 2s -m+ -m_\ 
V ™ 2s / 



m- ,m+ ,ni2s 



21 + p — 1 — m + — m 2s \ /7 — p + m_ — 1 



/2Z +p- 1 - m 2s - m+ - m_\ . 

<5(2to 2s + m_ + m + - 3p - 1) 



_ 2 (2l + p — m + — m 2s \ (I — p + ni 

^ * 

m- ,m+ ,ni2s 

f 2l + p — m 2s — m + — m 



(5(2m 2s + 77i_ + m + — 3p — 2) > . 
We reorganise the binomial coefficients as in (A. 3), to get 

^o(0-E{2 2 ^E( 2 ^ 2?, TO _ 1+m - 

p— o I. m— 



+ 22(i _ p) _ 2 ^^-2p-2 + m_j x (A19) 



I — p + m_ — 1\ (21 — 2p — 2 + m 2s + m 
2Z - 2p - 2 + m 



<5(2m 2s + m_ + m + — 3p — 1) 



+ 2 2(i-p)-2^ 



^ (l-p + m\ (2l+p-m 2s -m + -m\ 1 
> i 2m 2s + m_ +m + - 3p- 2) 

^ V m + / V TO 2s / I 

use the generating functions (A. 5) to perform the sum over m + and m 2s by picking 
out the appropriate residue, collapse the Kroenecker S and then use (A. 7) to get 

1 <T r^/(l + a:)(l-2x) 2 \f>l / 4 l \ 



,(l + x)Hl-2^, 
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{i 



(l-2ar) (l + x)(l 



Ax 2 



4x 3 



_2z) j 



r) 



27ri J V3x 3i (1 + - 2x) 2/ 7 (1 - 3a;) 
x |4a; 2 + x(l - 2x) + (1 + x)(l - 2a;)} 



,3/-l 



jM-1 



(A20) 



• M = l(mod 3) 

Here also, N(m+, ra_, m 2s ) is given by (3.18). We set M = 3/ + 1 and m P = 3p,3p- 2, 
3p - 1 for m P = 0, 1 and 2 respectively. Thus we set up (3.22) as 

tf(0 = E E 2 2 (^)^ +p -J l2s - m+ )( Z - p + 1 + +m -)x (A21) 



'2Z + p — m 2s — m + — m 
m 2s 



5(2rri2s + m- + m + — 3p) 



+ E 2 2( '~ p)+2 

p=l rrt- ,m+ ,ni2 s 



(2l+p—\ — m + - m 2s \ fl — p + m- + 1 



21 + p — 1 — m2s — rn + — m_ 



<5(2to2 S + m_ + m + — 3p + 2) 



/ 2? + p — 1 — m 2s — m + — m 



(21 + p — 1 - to + - to 2s \ (l — p + 1 



<5(2to 2s + m_ + 77i + — 3p + 1). 

to 2s 

Going through the same steps as in the case Q = 1, M = 0, we get 

S 1 1 (/) = 3 3 '. (A22) 
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• M = 2(mod 3) 

Here, N(m + ,m-,m 2s ) is given by (3.17). For this sector the sum to be evaluated 
is identical to the one in Q = 0, M = 2. Therefore, for M = 31 + 2, 

Si{l) = 3 3m = 3 M " 1 . (A.23) 

Q = 2 

The sums are done in a fashion completely analogous to those in Q = 1. 
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Table 1. The allowed values of A + XB, P a , Pb, m E , mp, 3m E + 2mp + P a + Pb, 
for the superintegrable 3-state chiral Potts model. 



M ee (mod 3) 


Q 


A + XB 


Pa 


Pb 


m E 




3m E + 2m P + P a + P b 


M~P a -P b -m P 














ee (mod 2) 


ee (mod 3) 


2M 


ee (mod 3) 


i 


-3- A 
-4A 
+3- A 
+2A 


2 
2 







1 
1 


ee 1 (mod 2) 
ee (mod 2) 
ee 1 (mod 2) 
ee (mod 2) 


ee (mod 3) 
ee l (mod 3) 
ee (mod 3) 
ee 2 (mod 3) 


2M - 1 
2M-2 
2M-2 
2M- 1 


ee 1 (mod 3) 
ee (mod 3) 
ee 2 (mod 3) 
ee (mod 3) 


2 


+3 + A 
+4A 
-3 + A 
-2A 




1 
1 


2 
2 




ee 1 (mod 2) 
ee (mod 2) 
ee 1 (mod 2) 
ee (mod 2) 


ee (mod 3) 
ee l (mod 3) 
ee (mod 3) 
ee 2 (mod 3) 


2M - 1 
2M-2 
2M-2 
2M- 1 


ee 1 (mod 3) 
ee (mod 3) 
ee 2 (mod 3) 
ee (mod 3) 



M = 1 (mod 3) 


Q 


A + XB 


Pa 


Pb 


m E 




3m E + 2m P + P a + P b 


M -P a -P b -mp 





-2 - 2A 
+1 + A 


2 





1 


= (mod 2) 
ee l (mod 2) 


ee (mod 3) 
ee 2 (mod 3) 


2M 
2M 


= 2 (mod 3) 
ee 1 (mod 3) 


l 


-2 

+1 - 3A 
+4 

+1 + 3A 




2 
2 




2 
2 


ee (mod 2) 
ee 1 (mod 2) 
ee (mod 2) 
ee 1 (mod 2) 


ee (mod 3) 
ee 1 (mod 3) 
ee 2 (mod 3) 
ee 1 (mod 3) 


2M- 1 
2M- 2 
2M~2 
2M - 1 


ee (mod 3) 
ee 2 (mod 3) 
ee (mod 3) 
ee 1 (mod 3) 


2 


-2 + 2A 
+1-A 










ee (mod 2) 
ee 1 (mod 2) 


ee (mod 3) 
ee 2 (mod 3) 


2M- 2 
2M- 1 


ee 1 (mod 3) 
ee 2 (mod 3) 




M ee 2 (mod 3) 


Q 


A + XB 


Pa 


-P& 


m E 


mp 


3m E + 2m P + P a + P b 


M -Pa-P,,- m P 





-1 - 1A 
+2 + 2A 


1 





2 


ee 1 (mod 2) 
ee (mod 2) 


= (mod 3) 
ee 1 (mod 3) 


2M 
2M 


ee 1 (mod 3) 
ee 2 (mod 3) 


l 


-1 + 1A 
+2 - 2A 










ee 1 (mod 2) 
ee (mod 2) 


ee (mod 3) 
ee l (mod 3) 


2M- 1 
2M-2 


ee 2 (mod 3) 
ee 1 (mod 3) 


2 


-4 

+1 - 3A 
+2 

-1 + 3A 


2 
2 






1 
1 


ee (mod 2) 
ee l (mod 2) 
ee (mod 2) 
ee 1 (mod 2) 


ee (mod 3) 
ee 2 (mod 3) 
ee 1 (mod 3) 
ee 2 (mod 3) 


2M-2 
2M- 1 
2M- 1 
2M-2 


ee (mod 3) 
ee 1 (mod 3) 
ee (mod 3) 
ee 2 (mod 3) 
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Table 2. The number of allowed sets of Vk (denoted by TV(sets) ) which 
solve (2.22) and the corresponding number of eigenvalues, TV (sets) x 2" lE , 
for each of the allowed sets of P a and Pb given in table 1 . 













M 


= 3 


Q 


Pa 


Pb 


TOE 


top 


TV(sets) 


iV (eigenvalues) = 2 E x iV(sets) 











2 





1 


4 











3 


5 


5 


total 


9 


1 


2 





1 





1 


2 




2 








1 


1 


1 







1 


1 





1 


2 







1 





2 


4 


4 


total 


9 


2 





2 


1 





1 


2 







2 





1 


1 


1 




1 





1 





1 


2 




1 








2 


4 


4 


total 


9 













M 


= 4 




Pa 


Pb 


TOE 


top 


TV (sets) 


iV(eigenvalues) = 2 mE x iV(sets) 





2 





2 





1 


4 











3 


5 


5 







1 


1 


2 


9 


18 


total 


27 


1 


1 





2 





1 


4 











3 


8 


8 




1 





1 


1 


3 


6 







2 





2 


3 


3 







2 


1 


1 


3 


6 


total 


27 


2 








2 





1 


4 











3 


5 


5 










1 


2 


9 


18 


total 


27 
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Table 2 (cont'd.) 













M 


= 5 


Q 


Pa 


Pb 


m E 


mp 


TV(sets) 


/^(eigenvalues) = 2™ E x AT(sets) 





1 





3 





1 


8 








1 


3 


23 


46 







2 


2 


1 


4 


16 











4 


11 


11 


total 


81 


1 








3 





1 


8 








1 


3 


23 


46 










2 


1 


4 


16 











4 


11 


11 


total 


81 


2 


2 





2 





1 


4 











3 


5 


5 




2 





1 


2 


9 


18 







1 


2 


1 


5 


20 











4 


16 


16 







1 


1 


2 


9 


18 


total 


81 













M 


= 6 


Q 


Pa 


Pb 


m E 


mp 


TV(sets) 


iV(eigenvalues) = 2™ E x iV(sets) 











4 





1 


16 








2 


3 


46 


184 











6 


43 


43 


total 


243 


1 


2 





3 





1 


8 








1 


3 


23 


46 




2 





2 


1 


4 


16 











4 


11 


11 







1 


3 





1 


8 








1 


3 


23 


46 







1 


2 


2 


19 


76 











5 


32 


32 


total 


243 


2 





2 


3 





1 


8 








1 


3 


23 


46 







2 


2 


1 


4 


16 











4 


11 


11 




1 





3 





1 


8 








1 


3 


23 


46 




1 





2 


2 


19 


76 











5 


32 


32 


total 


243 
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Table 3. Comparison for a six-site chain of the number of states of the 5=1, XXZ model with 
the number of allowed solutions of the superintegrable chiral Potts equation. The correspondence of 
momenta is that for mp = 0(mod3), Ps=i = Pcp and for mp = 1, 2(mod3), Ps=i — Pcp — 2n/3. 



chiral Potts 


S=l XXZ 


Ps=i 


Ps=i 


mp 


Q 


Pa 


Pb 





^ 6 


1 4-7T 

^ 6 


7T 


total 


s z 





_|_27T 

^ 6 


1 4-7T 

^ 6 


7T 


total 














1 











1 


6 


1 











1 


1 


i 


2 





1 


1 





1 


4 


5 


1 


1 


1 


1 


6 


2 


i 





1 


4 


3 


3 


3 


19 


4 


4 


3 


4 


3 


21 


3 











5 


8 


8 


9 


46 


3 


9 


8 


8 


9 


50 


4 


1 


2 





4 





3 


1 


11 


2 


16 


14 


16 


14 


90 


5 


1 





1 


6 


5 


5 


6 


32 


1 


21 


21 


21 


21 


126 


6 











10 


6 


8 


5 


43 





26 


21 


24 


23 


141 
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Table 4. The eigenvalues of the S — 1, 7 = 7r/3 XXZ periodic spin chain of six sites. The 
eigenvalues for negative values of S z are obtained using the symmetry of S z — > — S z . The underlined 
eigenvalues are obtained by using the solutions of (2.29) allowed for the superintegrable chiral Potts 
model in the energy formula (4.4). 



s z 


P 


ENERGY EIGENVALUES 


6 





18 


5 





15 




±27r/6 


16 




±4tt/6 


18 




7T 


19 


4 





12.4366930949376, 15.42532835575841, 19.10481891440097, 21.03315963490304 




±2tt/6 


13.94591578461395, 17.11126415759022, 19.94282005779584 




±47r/6 


14.19806226419517, 15.55495813208738, 17.24697960371747, 18 




7T 


16,16,16 


Q 
O 


n 


10.50604079256506, 13.89008373582526, 16, 16, 16.60387547160968, 

18 18 18 18 

lo, lo, lo, lu 




±2tt/6 


12.15498460202909, 14.82874774570659, 15.16848243075719, 15.6010235988223, 
16.39000960599962, 17.91683389294189. 18.09653253799185. 19.84338558575142 




±4tt/6 


13.13859446481224, 13.29978403077409, 14.30123932226858, 15.52206260527575, 
16.49825145060048, 16.68317951010723, 20.07618217185878, 20.48070644430282 




7T 


12.6791155005525. 12.67911550055252, 13.60074370344936. 14.64806079844653. 
14.64806079844655, 17.47824655651311, 18.67282370100094, 18.67282370100096, 
21.92100974003754 


2 





9.22646217887607, 12.43669309493758, 14.26794919243112, 14.26794919243113, 
14.49943450055038. 15, 15. 15, 15. 15.4253283557584, 16.84789556123169. 
17.73205080756888, 17.7320508075689, 19.10481891440096, 
21.03315963490305,23.42620775934184 




±2tt/6 


10.76229969320369, 13.12119617997307, 13.46167018036698, 13.94591578461395, 
14.82774302161002, 15.166048496364, 16, 16, 16.78684626768898, 
16.80912962588479, 17.11126415759021, 18.81398139161249, 
19.94282005779584,21.25108514329594 




±4tt/6 


11.50761996219082, 11.65982702674793, 12.20871215252208, 13.93039354102619, 
14.15544528228914, 14.19806226419516. 15.37778436506803. 15.55495813208736. 
16.79128784747792, 17.24697960371747, 18, 18, 18, 18.96238860818403, 
19.02049822873401,21.38604298575983 




7T 


12, 12.55051025721682, 12.55051025721683, 13.62771867673099, 
13.62771867673099, 16, 16, 16, 17.44948974278317, 17.44948974278318, 
19, 19, 19.37228132326903, 19.37228132326903 
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Table 4 (cont'd.) 



s z 


p 


ENERGY EIGENVALUES 


1 





8.5214914003045, 11.35042411530895, 12.43669309493757, 12.43669309493761, 
13.37655732137895, 14.26794919243112, 14.26794919243113, 15, 15, 15, 
15.42532835575839, 15.42532835575842, 17, 17, 17.73205080756889, 
17.73205080756892, 18.75152716300766, 19.10481891440092, 
19.10481891440093, 21.03315963490303, 21.03315963490304 




±2tt/6 


9.8720197837459, 10.76229969320372, 12.36017677219731, 13.12119617997308, 
13.46167018036698, 13.94591578461394, 13.94591578461395, 14.82774302161003, 
15.16604849636401, 15.52216410833847, 16, 16.08917517655664, 
16.78684626768899, 16.80912962588481, 17.11126415759022, 17.11126415759023, 
18.8139813916125, 19.15646415916171, 19.94282005779584, 19.94282005779586, 
21.25108514329595 




±4tt/6 


10.7817578014951, 11.50761996219081, 12.20871215252207, 12.77589846221546, 
13.93039354102618, 14.15544528228914, 14.19806226419515, 14.19806226419517, 
14.77231645874024, 15.55495813208737, 15.55495813208738, 16.79128784747792, 
17.24697960371746, 17.24697960371748. 17.38267160768822, 18, 18. 
18, 19.02049822873401, 21.28735566986096, 21.38604298575983 




7T 


11.13933320993551, 12.55051025721681, 12.55051025721682, 13, 13 
13.62771867673097. 13.62771867673099. 14.59566258403823. 
16, 16, 16, 16, 16, 16, 17.44948974278319, 17.44948974278321, 
18.51366656949581, 19, 19.372281323269, 19.37228132326903, 
23.75133763653045 








8.29991623268163. 10.50604079256506. 10.50604079256507. 10.94098877100413. 
13.08651430707976, 13.89008373582524, 13.89008373582526, 14, 14, 
14.76257547310441, 15.61946171799351, 16, 16, 16, 16, 
16.60387547160966, 16.60387547160969, 18, 18, 18, 18, 18, 18, 
19.19636205699761, 21.74114623721635, 24.35303520392256 




±2tt/6 


9.56483858138269, 12, 12.15498460202909, 12.1549846020291, 
12.76393202250021, 14.69952946524302, 14.82874774570659, 14.82874774570659, 
15.1684824307572, 15.16848243075721, 15.60102359882231, 15.60102359882231, 
16.39000960599962, 16.39000960599963, 17.23606797749979, 17.91683389294189, 
17.91683389294189, 18.09653253799186, 18.09653253799186, 19.84338558575142, 
19.84338558575144,21.73563195337429 




±4tt/6 


10.53589838486224, 11.22809490982445, 12.35288805589831, 13.13859446481223, 
13.13859446481224, 13.29978403077409, 13.29978403077409, 14.24071978009154, 
14.30123932226858, 14.30123932226858, 15.52206260527575, 15.52206260527575, 
16.49825145060049, 16.49825145060049, 16.55312816507534, 16.68317951010723, 
16.68317951010724, 17.46410161513774, 18.21877692510021, 19.40639216401014, 
20.07618217185877, 20.07618217185878, 20.48070644430281, 20.48070644430282 




7T 


10.87689437438233, 12.6791155005525, 12.67911550055251, 12.67911550055251, 
12.67911550055252, 13.60074370344935, 13.60074370344937, 14, 
14.64806079844652, 14.64806079844654, 14.64806079844654, 14.64806079844655, 
16, 16, 17.4782465565131, 17.47824655651312, 18.67282370100094, 
18.67282370100094, 18.67282370100095, 18.67282370100095, 
19.12310562561766, 21.92100974003751, 21.92100974003753 



39 



